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Abstract. 

The last 10-20 years has seen a profound shift in views of how the 
Galaxy's halo formed. The idea of a monolithic early collapse of a single 
system (Eggen, Lynden-Bell and Sandage 1962) has been challenged by 
observations at high redshift and by cosmological models of structure 
formation. These findings imply that we should see clear evidence of 
hierarchical formation processes in nearby galaxies. Recent studies of 
our Galaxy, made possible by large-scale CCD surveys such as the Sloan 
Digital Sky Survey (SDSS), have begun to reveal tantalizing evidence 
of substructure in the outer halo. We review evidence for tidal streams 
associated with known Milky Way satellites and for star streams whose 
progenitors are still unknown. This includes results from the SDSS and 
our own ongoing pencil-beam halo survey, the Spaghetti survey. 



1. Introduction 

A revolution is underway in our understanding of how galaxies form: The idea 
of a monolithic collapse of a single system at earliest times (Eggen et al. 1962) 
has been challenged by observations at high redshift and by increasingly so- 
phisticated models of the evolution of structures within the basic framework of 
the Hot Big Bang model (for example, Pearce et al. 1999, Steinmctz & Navarro 
1999). These findings imply that we should see clear evidence of hierarchical for- 
mation processes within nearby galaxies. Indeed, some of the earliest work that 
suggested complex and possibly hierarchical evolutionary histories for galaxies 
such as ours came from studies of the stellar populations of the outer halo of the 
Milky Way (Searle and Zinn 1978). Studies of nearby and high-redshift systems 
represent complementary approaches. 

We can learn much about accretion processes in galaxy formation using 
the Milky Way because we can view our halo in 3 dimensions. By adding the 
kinematics of halo stars to their 3-D positions, moving groups arc given enhanced 
contrast (see, for example, Harding et al. 2001). This sort of study is impossible 
for any other spiral galaxy. 

Recent studies of our Galaxy have begun to reveal tantalizing evidence of 
possible substructure in the halo, both locally (Majewski 1992, Cote et al. 1993, 
Arnold and Gilmore 1992, Helmi et al. 1999) and at larger distances (Ivezic et 
al. 2000, Yanny et al. 2000). And of course, the Sgr dwarf galaxy (Ibata et al. 
1995) is injecting globulars and field stars that formed in a completely separate 
entity from the Milky Way into its halo. The question has now shifted from 'Did 
accretion play any role in the formation of the Galaxy?' to 'How much of the 
Galaxy formed from accretion of hierarchical fragments?' We will show below 
that the data that is needed to answer the second question definitively also has 
the potential to teach us about the dark matter halo of the Galaxy and may 
directly constrain cosmological formation models such as that of Klypin et al. 
(1999), Moore et al. (1999). 
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2. The "Spaghetti" Survey 

Several years ago we began a search designed to find halo streams if they are 
present. Since then, the Sloan Digital Sky Survey (SDSS) team has found one 
stream in their first few hundred degrees of equatorial test data (Ivezic et al. 
2000, Yanny et al. 2000). Our group has found evidence that has strengthened 
the case that this substructure is associated with the Sgr dwarf using velocities 
of halo giants (Dohm-Palmer et al.2001). Our survey is unique in using a set of 
halo tracers (giants) that can be observed spectroscopically out to distances of 
hundreds of kpc where it is likely that only dwarf galaxies are found. 

Our photometric survey consists of several hundred pencil beams at high 
Galactic latitude, which will comprise 100 square degrees of data in the Wash- 
ington system. Our goal is to use this survey, the spectroscopic followup, and 
our extensive modelling, to quantify the amount of halo substructure associated 
with completely destroyed satellites, with partially destroyed satellites like the 
Sgr dwarf, and with intact objects such as to Cen. More details of the survey are 
available in Morrison et al. (2000), Dohm-Palmer et al.(2000), Morrison et al. 
(2001) and at http:/ /www. astro. lsa.umich.edu/users/rdpalmer/spag/] . In this 



paper we focus on the part of the survey which identifies distant halo giants. 

Local studies of halo stars passing through the solar neighborhood are lim- 
ited in numbers of stars, do not well sample all regions of the halo, and often 
have biased kinematics. Chiba and Yoshii(1998) have shown that the 250 red 
giants and RR Lyrae stars with space motions determined from Hipparcos oc- 
cupy a region extending to at most 35 kpc from the Galactic center. Carney et 
al. (1996) derive the apogalacticon distances of approximately 1200 local stars. 
There are only 23 stars in this sample with calculated apogalacticon distances 
larger than 30 kpc. These particular stars, by definition, are on extremely elon- 
gated orbits. There are vast parts of the Milky Way, and vast parts of phase 
space, therefore left unsampled. Figure 1 illustrates this clearly: halo globular 
clusters are plotted with open symbols, and it can be seen that they extend 
out to almost 100 kpc from the Galactic center. Known field stars (taken from 
the compilation of Beers &; Sommer-Larsen 1995) are shown with closed circles. 
Most halo field stars known are within a few kpc of the Sun. 

Figure 2 shows the spectrum of one of the most distant giants we have 
identified so far, star 1304.49b+60.51, which lies at a distance of ~80 kpc from 
the Sun. We expect to identify several hundred giants from the extreme outer 
halo when our 100 square degree survey is complete: an order of magnitude 
more than are currently known. 



3. Searching for Substructure in the Outer Halo 

Halo red giants are useful probes because of their high luminosity, which makes 
spectroscopic followup much easier. However, they are rare: we average ~1 giant 
in each 0.25 square degree field. This means that we cannot use many of the 
techniques which were developed to search for substructure in the larger local 
samples. These techniques are often one- and two-dimensional variations on the 
use of velocity histograms, whether to search for a spike (Cote et al. 1993), 
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Figure 1. Halo field stars known in 1995, from the compilation of 
Beers &; Sommer-Larsen (1995), plotted with closed circles. Compari- 
son with the distribution of globular clusters (open circles) shows that 
even though halo field stars are likely to stretch out to ~100 kpc from 
the Sun, very few are known in the outer halo. 
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Figure 2. Spectrum of one of the most distant halo giants wc have 
identified, star 1304.49b-|-60.51, which has a metallicity of [Fc/H]~-1.9 
and a distance of ~80 kpc from the Sun. Shown for comparison is the 
spectrum of a giant of similar color in the most metal-poor globular 
cluster known, NGC 5053 D. The features that we use for luminosity 
discrimination, Cal 4227 and the Mgfe/MgH feature, are marked. 
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unusual outliers (Majewski 1992) or simply a non-gaussian velocity distribution 
(Harding et al. 2001). 

We have used the extensive simulations of Harding et al. (2001), re-normalized 
so that one particle corresponds to one halo giant, to help us identify and test 
techniques to search for substructure. Our method uses the fact that that there 
are strong correlations between position and velocity along a disrupting stream, 
illustrated in Figure 3. Our sparse sampling reduces the numbers of stars de- 
tected in a given stream, but preserves the position/ velocity correlations, as can 
be seen in Figure 3. 

Since we can measure the distances to our giants with reasonable accuracy, 
we are able to derive all three spatial coordinates, but because they are too 
distant for proper motions until space missions such as GAIA and SIM are oper- 
ating, we only have one of the three velocity coordinates, the line-of-sight radial 
velocity Vios- Wc define a scaled metric in these four dimensions to measure the 
correlations between distance and velocity: two stars with similar distances and 
velocities will have a very small 4-D "distance" apart. 

For two stars with galactic latitude, galactic longitude, distance and velocity 
h,bi,di, Vi and l2,b2,d2,V2, we write the metric as: 

1 /2 

Dibdv = (wi{h - hf + W2{bi - bif + W3{di - dif + w^iV^ - Vs)^) 

We choose the scalings Wi for the 1, b, distance and velocity terms to account 
for the different measurement accuracy of each quantity, and for the fact that 
streams can spread spatially but have a tight relation between position and 
velocity at a given location. We know 1 and b very accurately, velocity to ~10% 
(20 km/s) and distance to ~25%, so two values that differ by less than their 
errors will be given less weight than values that differ significantly. 

We use two different simulated datasets to test this technique, with sample 
sizes equal to our expected number of giants when we have completed our survey. 
The first contains 500 points from an accreted halo simulation constructed with 
23 disrupting tidal streams from the library of Harding et al. (2001). These 
points were chosen from a dataset containing 320,000 points in 23 streams in 
a way that simulates our real observations. Wc used our real field centers for 
the 50 square degrees where observations have been taken, and planned field 
centers plus reflections about the galactic plane to make 400 pencil-beam field 
centers at high galactic latitude. If the stream particle was in the 0.25 square 
degrees surrounding a field center, we applied observational errors of 20 km/s 
for velocity and 25% for distance. The second dataset contained points from 
a smooth halo constructed with an isotropic velocity ellipsoid, as described in 
Harding et al. (2001), chosen to simulate our real observations in the same way. 

We calculate the value of the scaled 4-d metric described above for every 
pair of points in the smooth and accreted halo simulations. As expected, the 
accreted halo shows significantly more points with small values of the metric, as 
can be seen in Figure 4. We are working with Jiayang Sun (CWRU statistics) 
to develop sensitive statistical tests for this accretion signature. 
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Figure 3. Disruption of a stream in position and velocity space. The 
top panel plots the X and Y galactic coordinates of points along several 
wraps of the stream. The middle panel shows the line-of-sight velocity 
plotted against X, illustrating the velocity gradients along the stream. 
The bottom panel shows the particles that would be observed in a 100 
square degree survey such as ours, with realistic errors added in both 
distance and velocity. Although the picture is less clean, correlations 
between position and velocity still persist. 
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Figure 4. Comparison of 4-d metric values for a smooth halo (filled 
histogram) and for an accreted halo (open histogram). (Cross-hatched 
histogram shows particles matched with particles from another stream) . 
It can be seen that the accreted halo has significantly more points with 
small 4-d distances, as expected. 



4. First Results 

Here we show some preliminary results for our sample of 58 halo giants with 
spectroscopic confirmation. The distances are calculated using the photometric 
metallicities only at this point, although we will soon have spectroscopic metal- 
licities available, which will lead to more accurate distance measures. Figure 5 
shows the histogram of metric values for this sample (calculated as described 
above) and the histogram of metric values for a smooth halo sample with the 
same size as our real halo giant sample, chosen from the simulated smooth halo 
described above. While the sample size is small, we are already seeing some of 
the signatures of an accreted halo. 

Some of the smallest values of the 4-d metric in our giant sample come 
from a group of stars which are located in the same region of the sky as the 
SDSS overdensity (Ivezic et al. 2000, Yanny et al. 2000). This gives us the 
opportunity to test suggestions that this stream is associated with the Sgr dwarf 
by comparing the velocities of giant stars at similar distances to the SDSS BHB 
and RR Lyrae stars to the predictions of models of the disruption of Sgr. We 
have chosen to use the models of Helmi & White (2001), which, unlike earlier 
models of its disruption, can explain the fact that the galaxy was not completely 
disrupted on its first few passages without resort to unusual scenarios such as a 
very stiff and extended dark matter halo (Ibata & Lewis 1998) or a convenient 
recent collision with the Large Magellanic Cloud (Zhao 1998). 

Figure 6 shows the predictions of the purely stellar model of Helmi Sz White 
(2001) for particles close to the minor axis of the Galaxy, compared with the 
actual positions and velocities of 21 confirmed metal-poor giants in this direction. 
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Figure 5. Comparison of 4-d metric values for a smooth halo (filled 
histogram) and for our preliminary halo giant sample (open histogram). 
The smooth halo was constructed to have the same sample size as our 
real sample. Although the sample size is small, we are already seeing 
indications of some clustering in the 4-d velocity/distance space. 



It can be seen that all of the halo giants with distances greater than 30 kpc agree 
well with the predictions of the Sgr disruption model, confirming the suggestions 
of Ivczic ct al. (2000) and Ibata ct al. (2001a) that the SDSS stream is associated 
with Sgr. The right-hand panel of the plot shows, for comparison, one simulation 
taken from a smooth halo: it can be seen that the arrangement of distances and 
velocities that wc observe is very unlikely to have occurred by chance. More 
details of these observations can be found in Dohm-Palmer et al. (2001). 

The wrap at distances near 50 kpc is the one detected by the SDSS. It is 
interesting to note that we may also have detected stars from the wraps near 
20 kpc (although the large velocity range predicted for these wraps makes this 
harder to test) and from the wrap at 80 kpc, which contains debris lost from the 
Sgr dwarf on one of its earliest passages. Such debris, if confirmed with a larger 
sample, will give us the opportunity to measure the M/L ratio of the progenitor 
to the Sgr dwarf, as the outer regions arc the first to be lost to tidal disruption, 
and the presence of a massive dark matter halo in the outer regions will help 
retain more stars than in a galaxy with little or no dark matter (Helmi 2000). 

Our sample of giants from 100 square degrees of imaging will be large enough 
to assign membership to different streams in the outer halo. This will allow 
their use in ways that were almost inconceivable ten years ago. For example, 
we can learn about the structure of the progenitor satellites: the red giants that 
are Sgr stream members have a mean metallicity of [Fe/H]=-1.5, significantly 
more metal-poor than the main body of Sgr (Layden and Sarajedini 2000). 
This suggests that the Sgr progenitor had a radial abundance gradient. As 
mentioned above, identification of stars lost on early passages of Sgr will allow 
us to constrain the mass-to-light ratio of its progenitor. The spread in longitude 
of the Sgr debris will constrain the flattening of the Galaxy's dark matter halo 
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Figure 6. (Left panel) Comparison of the purely stellar model of the 
Sgr dwarf and its disruption from Helmi & White (2001) (closed black 
symbols) with the 21 halo giants we have found in this direction (grey 
symbols with error bars) . The right panel shows the comparison of the 
model predictions with a typical simulation of a smooth halo. It is very 
unlikely that our data came from a smooth halo. 



(Helmi 2000, Ibata et al. 2001b). We can even use the stream velocity dispersion 
to constrain the lumpiness of the Galaxy's dark halo, which will allow us to 
directly test the predictions of simulations of structure formation (eg Klypin et 
al. 1999, Moore et al. 1999). 
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